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^ H We show for a first time ever a prototype of a fully exclusive QCD NLO parton shower for 

Q-i the initial state (albeit for a limited set of diagrams). It is based on the rigorous theorems of 

the collinear factorisation, however the standard DGLAP evolution kernels are replaced by their 
!— 1 exclusive versions. Contrary to the standard DGLAP approach, the constructed parton shower 

provides fully exclusive events, i.e. four-momenta of all emitted partons. At the inclusive level it 

is identical to the MS scheme of the DGLAP evolution. 
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With the approach of the first data from the LHC experiments the issue of the precision of the 
Monte Carlo (MC) simulations of the perturbative QCD becomes more and more important. In the 
domain of the semi-inclusive, analytical or semi-analytical calculations the results are spectacular: 
the NLO accuracy is a standard and numerous quantities are calculated to much higher precision, 
just to mention the DGLAP kernels known to NNLO level [1]. On the contrary, all the currently 
available MC Parton Shower (PS) codes are of the LO or, at most, of the improved LO type (as 
far as the differential distributions are concerned, the overall, inclusive, normalization can easily 
be corrected to higher order precision), based on the methods developed in early 1980-s, cf. for 
example [2, 3, 4]. 

It is therefore justified to ask whether it is possible to develop a new, more precise scheme 
for the QCD MC parton shower? In this presentation we will demonstrate that yes, it is feasible. 
We will describe a working prototype of such a genuine NLO parton shower for the initial state 
QCD which we have developed recently within the KRKmc project. Note that there are some other 
attempts to construct parton shower schemes beyond leading order as well, cf. eg. [5, 6, 7]. 

The KRKmc parton shower fulfills the following list of requirements: 

(A) it is based firmly on Feynman Diagrams (Matrix Element) and Lorentz-invariant phase space 
(LIPS); 

(B) it is based rigorously on the collinear factorization [8, 9]; 

(C) it implements exactly the NLO MS DGLAP evolution [10]; 

(D) it implements fully unintegrated parton density functions (PDFs); with NLO evolution done by 
MC itself, using new, exclusive NLO kernels. 

The combination of the collinear factorization (item 1), which is the best proven factorization 
scheme, with the complete kinematical information of the generated partons (item 4), i.e. full ex- 
clusiveness of the shower is the essence of the novelty of our approach. Let us now highlight some 
of the important points in this new KRKmc scheme. For more details we refer to Refs. [11, 12, 13]. 

Step 1. Re-do the calculation of the NLO kernels, based on the framework of Curci, Furmanski 
and Petronzio [10] but in the exclusive way. Use various types of evolution time (virtuality, k±, 
rapidity . . . ). The standard DGLAP kernels do not depend on the choice of the evolution time. On 
the contrary, the exclusive ones depend. Use MS factorization scheme. In the Fig. 1 we show all the 
contributions to the non-singlet LO and NLO evolution kernels. In this work we restrict ourselves 
to the Cp- part (two "bremsstrahlung-like" graphs) labeled with the blue box. It is in a sense the 
most complicated part because it involves the subtraction of the soft counter-term (labeled as 1 — P 
in the plot). 

Step 2. Remove the dimensional regularization, go back to four-dimensions. Use instead a 
geometrical cut-off A for the collinear singularity. But, all the time stay within the MS scheme. 
This is possible because the collinear pole comes from a simple, factorisable integral of the form 




Step 3. Re-formulate the factorization formula, because the DGLAP equation mixes orders of 
perturbative expansion. The NLO DGLAP kernel, denoted as P, is in fact a sum of LO and NLO 
pieces. Therefore P k terms are a mixture of various orders: 

P = aP L0 + a 2 P NL0 P 2 = a 2 {P LO ) 2 + a\P LO P NLO + P NLO P LO ) + a\P NLO ) 2 . 
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Figure 1: Contributions to the NLO non-singlet evolution kernel. The dp part label-led with the blue box. 

Such a set-up, with partial, incomplete higher orders, may lead for example to the negative events 
in the MC and should be avoided. 

Step 4. Resign from the ordering in the evolution time in the underlying LO crude MC. Use 
the Bose-symmetric form instead: 

rtnm JL 1 rtma JL 

/ n*<«Wi Jn / (D 

The ordering is an approximate feature of the LO matrix element and therefore it is not strict at the 
NLO level and we have to explicitly sum over entire phase space. 

Step 5. Construct appropriate MC weight with the NLO exclusive kernel. It will be then 
applied on the top of the standard LO MC. 

Step 6. Resolve the mismatch of the lower limits of the "internal NLO phase space", i.e. of 
the internal degree of freedom integrated out in the construction of the inclusive kernel. In the 
analytical calculation of the MS kernel the lower limit of dti ntema i integration is set to 0, whereas 
in MC it is limited by some f m ; n for all partons. In this work we resolve this conflict by "artificial" 
lowering of the f m ; n limit, below the actual start of the evolution and by performing a LO pre- 
evolution in this extended region. This way we maintain exact agreement with the MS result. 

We have implemented the above scheme of the exclusive NLO evolution in the MC program 
KRKmc. In the Fig. 2 we compare its results with a standard NLO DGLAP MC evolution. Both 
evolutions are implemented as weights on the top of the same LO Markovian MC algorithm. The 
curves shown correspond to the contributions with one and two NLO "insertions". Evolution ranges 
from 10 GeV to 1 TeV. LO pre-evolution starts at lGeV from 5(1 — x) distribution. As one can 
see the agreement is very good, within the statistical errors. This way we demonstrated for the 
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Figure 2: Comparison of exclusive and standard DGLAP NLO evolutions 

first time ever that the QCD NLO parton shower can be constructed (although for a limited set of 
Feynman diagrams) ! 

We are currently in the process of adding the rest of graphs from Fig. 1 , omitted in this work. 
This way the non-singlet evolution will be completed. Once also the singlet evolution is added, we 
will be ready to construct the complete event generator for the Drell-Yan-type processes at LHC 
or DIS at HERA. 
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